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Determining parameters of an equivalent circuit representing a 
5 transmission section of an electrical network 



DESCRIPTION 

10 

Field of the Invention 

The invention relates to large-scale electric power transmission networks, and, 
more particularly, to a method, computer program and apparatus for 
15 determining parameters of an equivalent circuit representing a transmission 
section of an electrical network as described in the preamble of claim 1 , 9 and 
10. 



20 Background of the Invention 

Electric power transmission and distribution systems or networks comprise 
high-voltage tie lines for connecting geographically separated regions, and 
substations for transforming voltages and for switching connections between 
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lines. Power generation and load flow in a network with several substations is 
controlled by a central energy management system. Ah important issue in the 
control of a power generation and load flow is to keep the network stable, i.e. 
to avoid voltage collapses and swings. A method for assessing network 
stability, based on voltage margins, is described in EP 1 217 709. This patent 
application and articles cited therein describe a "Voltage Instability Predictor" 
(VIP) that measures currents and voltages locally in order to infer a proximity 
to voltage collapse. The concept of the VIP is based on an equivalent network 
as shown in Figure 3. One part of an electric power system is treated as a 
power source, another part as a load. The power source is represented by its 
Thevenin equivalent 1 with a source voltage and a Thevenin or source 

impedance Z th . An apparent load impedance Z, represents the load. Both the 
Thevenin impedance Z th and the apparent load impedance Zj are estimated 
from the current and voltage measurements i 2 ,v 2 by a VIP device. The relation 

of these impedances, expressed by a stability margin or power margin, 
indicates how close the power system or network is to collapsing. 

This approach uses measurements taken at a single load bus. A single set of 
measurements does not contain enough information to directly compute the 
parameters of the Thevenin equivalent E th9 Z th , however these can be 

estimated using for example using the least-squares method once at least two 
sets of measurements made at different times are available. For accurate 
estimation of these parameters, a sufficient change in the measurements must 
have occurred, e.g. due to load change. During this time the generation and 
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transmission network is assumed to remain constant, which introduces an 
estimation error. Furthermore, the estimation is noise-sensitive and introduces 
a time-delay. The practical applicability of the method is therefore limited. 

Other approaches are based on a complete set of system information, that is, 
all voltages and currents are measured, and the impedances of all lines and 
loads and generator capability diagrams are measured or estimated. Although 
these techniques are powerful tools capable of very accurately assessing the 
voltage stability of the network, the large number of measurements required 
makes them expensive and impractical to implement in practice. 

It is therefore desirable to provide improved and efficient means to determine 
a model of a section of an electrical network, in particular for a section that 
essentially serves for transmitting power from a generating section to a load 
section. 

Description of the Invention 

It is therefore an object of the invention to create a method, computer program 
and apparatus for determining parameters of an equivalent circuit representing 
a transmission section of an electrical network of the type mentioned initially, 
which overcomes the disadvantages mentioned above and allows a fast 
assessment of network stability. 
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These objects are achieved by a method, computer program and apparatus for 
determining parameters of an equivalent circuit representing a transmission 
section of an electrical network according to the claims 1 , 9 and 10. 

5 In the inventive method, the transmission section is representable as having at 
least two interfaces with other sections of the network, and the method 
comprises the steps of 

a) determining, for each of the interfaces, a voltage phasor at the interface 
and a phasor of a current flowing through the interface, the 

10 measurements at the different interfaces being made essentially 

simultaneously, and 

b) computing, from said voltage and current phasors, values of impedances 
constituting the equivalent circuit. 



15 This allows determining the equivalent circuit from a single set of essentially 
simultaneous measurements. There is no time delay incurred by having to 
measure a time series of phasor values and estimating impedances, as in the 
state of the art. 



20 In a preferred embodiment of the invention, at least one interface comprises at 
least two physical power lines, and the voltage phasor at the interface is 
determined as a weighted sum of the voltages at the power lines. This allows 
representing a transmission section that is connected to a load or generation 
section with more than just a single power line. 



25 
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In a further preferred embodiment of the invention, a current phasor 
representing a current through the interface is preferably computed from the 
voltage phasor at the interface and from a power flow through the power lines 
constituting the interface. 

The computer program for determining parameters of an equivalent circuit 
representing a transmission section of an electrical network according to the 
invention is loadable into an internal memory of a digital computer, and 
comprises computer program code means to make, when said computer 
program code means is loaded in the computer, the computer execute the 
method according to the invention. In a preferred embodiment of the 
invention, a computer program product comprises a computer readable 
medium, having the computer program code means recorded thereon. 

A data processing system or apparatus for determining parameters of an 
equivalent circuit representing a transmission section of an electrical network 
according to the invention comprises means for carrying out the steps of the 
method according to the invention. In a preferred embodiment of the 
invention, the data processing system is an apparatus comprising a data 
processor, a memory coupled to the processor and computer program code 
means stored in said memory, where said computer program code means, 
when executed by the processor, causes the method according to the invention 
to be executed. 
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Further preferred embodiments are evident from the dependent patent claims. 



Brief Description of the Drawings 

5 

The subject matter of the invention will be explained in more detail in the 
following text with reference to preferred exemplary embodiments that are 
illustrated in the attached drawings, in which: 



10 Figure 1 schematically shows a power network with a transmission corridor; 
Figure 2 shows an equivalent circuit of the power network with separately 

modelled generating and transmission sections; 
Figure 3 shows an equivalent circuit of the power network with a combined 
model for the generating and transmission sections; 
15 Figure 4 schematically shows a power network with a transmission section 
having three interfaces to other sections of the power network; 
Figure 5 shows an equivalent circuit of the transmission section of figure 4; 
Figure 6 schematically shows power network with a transmission corridor 
having more than one power line at interfaces to other sections of 
20 the power network; and 

Figure 7 shows time histories of impedances and Thevenin source voltage 
for an equivalent circuit of a collapsing network. 
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The reference symbols used in the drawings, and their meanings, are listed in 
summary form in the list of reference symbols. In principle, identical parts are 
provided with the same reference symbols in the figures. 



Detailed Description of Preferred Embodiments 

Figure 1 schematically shows a power network with a transmission section 4 
having a first interface 5 to a generating section 1 and a second interface 6 to 
a load section 2. The transmission section 4 has two interfaces 5,6 and is 
therefore also called a transmission corridor. The object of the invention is to 
determine parameters of an equivalent network or circuit that represents the 
electrical behaviour of the transmission corridor 4. In a preferred embodiment 
of the invention, the equivalent is a T-equivalent network comprising a 
transmission impedance Z T and a shunt impedance Z sh as shown in Figure 2, 
in another preferred embodiment of the invention, the equivalent is a n- 
equivalent. The generating section 1 is modelled separately from the 
transmission section 4 and represented by a Thevenin equivalent comprising 
E g » Z g • 

The transmission section 4 serves to transmit electric power from the 
generating section 1 to the load section 2 of the network. These sections are 
predominantly generating, transmitting and load sections. In reality it is 
possible for a generating section 1 to comprise loads, for the transmission 
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section 4 to comprise loads and generators, and for the load section 2 to 
comprise generators. However, effects of such loads and generators is 
assumed to be negligible with respect to overall power flow. These generators 
and loads are then automatically included in the shunt impedance of the 
5 equivalent circuit, and are thus assumed to exhibit approximately impedance 
characteristics. However, the equivalent network does not accurately represent 
a generator that is controlled to provide a constant voltage or power. In such 
cases, generated or consumed power of power generators and loads that are 
part of the transmission section 4 should be small in comparison to the power 
10 transmitted through the transmission section 4. 

, According to the invention, phasor measurements are made at both ends or 
interfaces of the transmission corridor 4. This allows splitting the estimation 
into two stages. First, the parameters of the equivalent of the transmission 

15 corridor are determined through direct calculation from measurements made 
at both interfaces and synchronised, i.e. at essentially the same time. Second, a 
Thevenin equivalent of the feeding generators is computed. Once the 
parameters of the transmission section 4 and generation section 1 are known, 
stability analysis can be carried out analytically and various stability indicators 

20 may be calculated using well-known methods. An advantage over the present 
state of the art is that parameters of the equivalent network are computed 
from essentially simultaneous phasor measurements, in contrast to VIP-type 
approaches that need a series of measurements over time and thus introduce a 
time delay for obtaining a network parameter estimate. 

25 
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Phasor data is determined with a phasor measurement unit (PMU) that resides, 
for example, at a feeder at the bay level of substations or at branching points 
along transmission lines. A voltage phasor represents, for example, a voltage 
of the feeder or line, while a current phasor represents current flowing through 
5 the feeder or line. 

The phasor data represents a phasor and may be a polar number, the absolute 
value of which corresponds to either the real magnitude or the RMS value of a 
quantity, and the phase argument to the phase angle at zero time. 

10 Alternatively, the phasor may be a complex number having real and imaginary 
parts or the phasor may use rectangular or exponential notation. Phasors may 
be used to represent quantities such as the voltage, current, power or energy 
associated with a phase conductor or an electronic circuit. By contrast, 
conventional sensing devices used in power networks generally measure only 

is scalar, average representations, such as the RMS (root mean square) value of a 
voltage, current etc. 



In some PMU applications, the phasor data is collected from phasor 
measurement units that are distributed over a large geographic area, i.e. over 

20 tens to hundreds of kilometres. For applications in which the phasor data from 
these disparate sources are analysed in conjunction, they must refer to a 
common phase reference. In other words, the different phasor measurement 
units have local clocks that are synchronised to within a given precision. Such a 
synchronisation of the phasor measurement units is preferably achieved with a 

25 known time distribution system, for example the global positioning (GPS) 
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system. In a typical implementation, the phasor data is determined at least 
every 200 or every 100 or preferably every 40 milliseconds, with a temporal 
resolution of preferably less than 1 millisecond. In a preferred embodiment of 
the invention, the temporal resolution is less than 10 microseconds, which 
corresponds to a phase error of 0.2 degrees. Each measurement is associated 
with a time stamp derived from the synchronised local clock. The phasor data 
therefore preferably comprises time stamp data. 

The parameters of a T-equivalent of the transmission corridor 4 are calculated, 

including any load or generation that may be present in the transmission 

corridor. Applying Ohm's and Kirchhoff s laws, with the measured phasors 
denoted by complex values v^ij and v 2 ,* 2 , the complex impedances Z T ,Z sh and 

Z L are calculated as follows: 

Zj* — 2 ~ z — — — 
i, - i 2 

(1) z sh-— — 

h ~ l 2 

The complex voltage E g and impedance of the equivalent voltage source Z g 

cannot be calculated in the same straightforward way, so one of them must be 
assumed to be known to avoid the time-delay of an estimation procedure 
similar to that in the VIP-type approaches. If the generators have voltage 
controllers and can be assumed not to have capability limits, E g can assumed 

to be constant and Z ? can then be calculated using the formula: 
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(2) Z g =~^j 

However, in most practical cases it is more realistic to assume that Z g is 

known since it typically comprises step-up transformers and short 
transmission lines leading to the first interface 5 to the transmission corridor. 
5 Since the lines are short, they are less prone to change. In a preferred 
embodiment of the invention, the equivalent complex voltage of the 
generators is therefore computed as: 

(3) E g =v l +Z g i l 

Once the parameters of the T- and Thevenin equivalent are calculated, a 
10 second Thevenin equivalent E lh ,Z th for the combined generation and 

transmission corridor according to figure 3 can be calculated: 



^ Z T 1 



Z.u = — + 



2 J-^ 1 



(4) Zsh ^- + Z g 

2 8 

- _ Z th +Z L 
E th =v 2 = 



Based on this second Thevenin equivalent E th ,Z th , stability analysis can be 

1 5 performed analytically in a known fashion, e.g. be detecting when the Thevenin 
Z th impedance becomes smaller than the load impedance Z L . In a preferred 

embodiment of the invention, for determining the stability of a power network, 
the steps described above for determining the Thevenin equivalent of the 
generating section 1 and transmission section 4 are performed, followed by a 
20 further step of detecting a sudden change in the Thevenin source voltage E th 
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as an indicator of network instability. Figure 7 shows an example of 
trajectories over time t of Thevenin impedance Z th in relation to load 
impedance Z L and of Thevenin source voltage E th . At the time t=30s a large 
load change causes oscillations in both load impedance Z L and in Thevenin 
source voltage E th , well before instability is reached. Shortly after t=60, Z L 
becomes smaller than Z th , which is a known indicator of instability. Shortly 
before t=60, a step change in E th indicates approaching instability. 

In a further preferred embodiment of the invention, the transmission section 4 
has three interfaces to other sections 1,2,3 of the network, as shown in figure 
4, where a further section 3 is connected to the transmission section 4 through 
a further interface. In this case, the transmission section 4 is preferably 
modelled by an equivalent circuit as shown in figure 5. The equivalent circuit 
comprises three line impedances Z XTr ,Z 2yT ,Z 3lT interconnecting the 
interfaces having voltages Vj,v 2 ,v 3 and three shunt impedances 
Zi shJ Z 2sh yZ 3sh connecting the interfaces to a common node. The values of the 

impedances are determined from voltage and current phasors at the interfaces 
by application of Kirchhoff s laws. In a similar way, the inventive idea can be 
generalised to networks with a larger number terminals. 

In a further preferred embodiment of the invention, at least one interface 
between the transmission section 4 and other sections comprises at least two 
physical or actual power lines, as is shown in figure 6. For a network as in 
figure 6, first main load and generation centres or sections must be identified. 
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fn this case, a distinct generating section 1 is found in the area above a first 
cut line labelled 5, which comprises three major generators and some shunt 
compensation but only a few minor loads. Between the first cut line 5 and a 
second cut line 6 lies a section with no generation equipment and only a few 
5 minor loads, which shall be considered as transmission section 4. This is the 
transmission corridor, whose stability is of interest. In the equivalencing 
procedure described above, the loads will be implicitly included in shunt 
impedances. 

10 Below the second cut line 6 is an area with predominantly load character. There 
are some minor generation units, but in cases where the voltage stability is 

j 

endangered, these generators would have exceeded their capability limits and 
thus no longer contribute to stabilization. It is therefore reasonable to include 
them in the shunt impedance modelling the load. 

15 

Network section interfaces or region boundaries are now defined by the two 
transfer cuts 5 and 6. According to a preferred embodiment of the invention, 

I two virtual buses are defined, one for each end of the transmission corridor 4. 

i 

Buses 6, 13 and 14 of the original system are grouped into a virtual bus 
| 20 corresponding to the first interface 5, and buses 24, 15 and 16 are grouped 
into a virtual bus corresponding to the second interface 6. The part of the 
system between cuts 1 and 2 is considered a virtual transmission corridor 
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The active and reactive power through the interfaces 5,6 is computed as the 
sum of active and reactive power transfer p,q on ail lines crossing the 

respective interface 5,6 as follows: 

PcutX ~ P\4_S + P6_8 + P6_9_l + />6_9_2 + Pl3_ll_l + />13_11_2 

, v Qcut\ - <?14_8 + #6„8 + #6_9_1 + #6_9_2 + 9l3_ll_l + #13_11_2 
5 (5) 

Pcutl ~ P24_8_l + />24_8_2 + Pl5_10_l + / 7 15_10_2 + P\6 _9 _1 + jPl6_9_2 
#c**2 = #24_8_1 + #24_8_2 + 3l5_10_l +#15_10_2 + 9l6_9__l + #16_9_2 

Since the voltages v x and v 2 are voltages of virtual buses, they cannot be 

directly measured and must be extracted from other measurements. In a 
preferred embodiment of the invention, they are computed as a weighted 
10 average of the voltage of the buses the interface comprises. The weightings 
are proportional to the active power transported by the lines originating at that 
bus through the relevant transfer cut. 
(6) 

_ P\4 8^14 + (/>6 8 + P6 9 9 2^6 + (P\3 11 1 + P\3 11 2>^3 

Vi = = = — — =— = = — 

Pl4^8 + P6_S + P6_9_1 + ^6_9_2 + Pl3_ll_l + Pl3_ll__2 

(P24 8 1 +P24 8 2)^24 + (Pl5 10 1 + P\5 10 2^15 + ( P\6 9 \ + P\6 _9 ^2^16 
y 2= ^ 

Pl4_S^l + / 7 24_8_2 + Pl5_10_l + Pl5_10_2 + Pl6_9_l + / ? 16_9__2 

15 

Other variants for computing v, and v 2 would be the mean value, or simply 
selecting one of the available voltage measurements from the area. With the 
power flow through the cuts and the voltages of the virtual buses known, in a 
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preferred embodiment of the invention, the currents at the interfaces 5,6 of 
the virtual transmission corridor 4 are computed as 



5 From these measurements, estimation of the parameters of the virtual 
transmission corridor can be computed using the procedure described above. 
The method of indirectly computing the resulting equivalent currents from the 
calculated powers makes it possible to combine current measurements from 
different voltage levels. If all the included lines have the same rated voltage, or 
10 all calculations are carried out in per-unit quantities, then the resulting 
currents can be computed by adding each of the currents of the lines that 
make up the corresponding interface. 

List of designations 
15 1 generating section 

2 load section 

3 further section 

4 transmission section 

5 first interface 

20 6 second interface 
7 third interface 



(7) 
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PATENT CLAIMS 

A method for determining parameters of an equivalent circuit 
representing a transmission section (4) of an electrical network, where 
the transmission section (4) is representable as having at least two 
interfaces (5,6,7) with other sections (1 ,2,3) of the network, 
characterised in that the method comprises the steps of 

a) determining, for each of the interfaces (5,6,7), a voltage phasor 
(vj,V2,v 3 ) at the interface (5,6,7) and a phasor of a current (iiJ 2 *h) 

flowing through the interface (5,6,7), the measurements at the 
different interfaces (5,6,7) being made essentially simultaneously, 
and 

b) computing, from said voltage (v x> v 2 ,v$) and current (if, 12*13) 
phasors, values of impedances (Z T9 Z sh ) constituting the equivalent 

circuit. 

Method according to claim 1 , wherein the transmission section (4) is a 
transmission corridor having exactly two interfaces (5,6,7) to other 
sections (1 ,2) of the network. 

Method according to claim 2, wherein a first interface (5) connects the 
transmission corridor (4) to a network section (1) consisting 
predominantly of power generators, and a second interface (6) connects 



r 
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the transmission corridor (4) to a network section (2) consisting 
predominantly of power consumers. 

4. Method according to claim 2, wherein the transmission network (4) is 
represented by one of a T-equivalent and a n-equivalent circuit. 

5. Method according to claim 3 or 4, comprising the further step of 
computing parameters (E thy Z th ) of a Thevenin equivalent of a network 

constituted by the transmission section (4) and by the network section 
(1) consisting predominantly of power generators. 

6. Method according to claim 1, wherein the transmission section (4) 
comprises three or more interfaces (5,6,7) to other sections of the 
network (1,2,3) and the equivalent circuit comprises line impedances 
(Z, 2r ,Z23 r ,Z 3]7 -) interconnecting the interfaces and shunt impedances 
(Z lsh ,Z 2sh ,Z 3sh ) connecting the interfaces (5,6) to a common node. 

7. Method according to one of the preceding claims, wherein at least one 
interface (5;6;7) comprises at least two physical power lines, and the 
voltage phasor (v,;v 2 ) at the interface (5;6;7) is determined as a 

weighted sum of the voltages at the power lines. 

8. Method according to claim 7, wherein a current phasor (I X J 2 ) 
representing a current through the interface (5;6;7) is computed from the 
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voltage phasor iy x \v 2 ) at the interface and a power flow 
(Pcuti^cutuPcuti^cuti^ through the power lines constituting the 
interface. 

5 9. Computer program for determining parameters of an equivalent circuit 
representing a transmission section of an electrical network which is 
loadable and executable on a data processing unit and which computer 
program, when being executed, performs the steps according to one of 
the claims 1 to 8. 

10 

10. Data processing system for determining parameters of an equivalent 
circuit representing a transmission section of an electrical network 
comprising means for carrying out the steps of the method according to 
any one of the claims 1 to 8. 
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ABSTRACT 

In a method, computer program and apparatus for determining parameters of 
an equivalent circuit representing a transmission section of an electrical 
5 network, the transmission section (4) is representable as having at least two 
interfaces (5,6,7) with other sections (1,2,3) of the network, and the method 
comprises the steps of 

a) determining, for each of the interfaces (5,6,7), a voltage phasor 
(v lt v 2 ,v 3 ) at the interface (5,6,7) and a phasor of a current (i,,i 2 ,i 3 ) 

io flowing through the interface (5,6,7), the measurements at the different 

interfaces (5,6,7) being made essentially simultaneously, and 

b) computing, from said voltage (vi,v 2 ,v 3 ) and current (r,,z 2 ,i 3 ) phasors, 
values of impedances (Z T ,Z sh ) constituting the equivalent circuit. 

15 This allows determining the equivalent circuit from a single set of essentially 
simultaneous measurements. 

(figure 1) 
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